The Ji — J2 antiferromagnet on the square lattice with Dzyaloshinskii-Moriya 
interaction : An exact diagonahzation study 
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We examine the influence of an anisotropic interaction term of Dzyaloshinskii-Moriya (DM) type 
on the groundstate ordering of the Ji — J2 spin-i-Heisenberg antiferromagnet on the square lattice. 
For the DM term we consider several symmetries corresponding to different crystal structures. 
For the pure Ji ~ J2 model there are strong indications for a quantum spin liquid in the region 
of 0.4~J2/Ji~0.65. We find that a DM interaction influences the breakdown of the conventional 
antiferromagnetic order by i) shifting the spin liquid region, ii) changing the isotropic character 
of the groundstate towards anisotropic correlations and iii) creating for certain symmetries a net 
ferromagnetic moment. 



I. INTRODUCTION 



The two-dimensional quantum spin Heisenberg antiferromagnet (AFM) has attracted a considerable interest in 
connection with the magnetic properties of the high-temperature superconductors in recent time [0| . The CuO planes 
\^ being responsible for the superconductivity show a strong in-plane exchange and only a small ofF-planc exchange [E|. 

, Therefore in the insulating phase the interacting Cu spins should be well described by a two-dimensional spin-^- 
Tj" ' Heisenberg antiferromagnet. There are several arguments 1^-^ that additional to the nearest neighbour exchange 
\ Ji a frustrating diagonal J2 bond is relevant. The groundstate properties of this so-called Ji — J2 model are widely 
discussed in the last years, mainly with respect to a possible breakdown of the magnetic long-range order (LRO) due 
to the combined influence of quantum fluctuations and frustration (see e.g. [Q- p^]). One finds evidence for a finite 
parameter region around J2/ Ji^ 0.5 where quantum fluctuations open a window of a spin liquid phase. Moreover, 
there are indications of enhanced exotic order parameters (spin-Peierls, chiral) ||8|, [l0| , ^ , p^Jl9| . 

However, there are also indications for additional anisotropic terms in the Hamiltonian which could explain the 



'"^ ■ experimentally observed weak ferromagnetism for instance in La2Cu04 21 1. In general, a small ferromagnetic 

^ ' moment in antiferromagnets may appear in materials with low crystal symmetry. This tilting of the spins can be 

O described by adding the so-called anisotropic Dzyaloshinskii-Moriya (DM) interaction to the isotropic Heisenberg 

J-} model. Already 1957 Dzyaloshinskii formulated a phenomenological theory of these facts ||2^. Three years later 

^ : Moriya developed the microscopic theory of the weak ferromagnetism [p3[ . The occuring additional interaction term 



\ in the Hamiltonian is proportional to the DM vector D. In some recent publications ||24|-|28| one has examined the 
5^ ' origin and the structure of D. In this paper we discuss the influence of DM terms on the groundstate properties of 
' the Ji — J2 model on the square lattice. 

II. ANISOTROPIC SPIN INTERACTION 

Let us start with a general Hamiltonian, describing quadratic spin-spin interaction: 

N 

H=Y, Si J„- S, (1) 
with Ji j being a 3x3 matrix of interaction, which can be written as 



J ^\{SpJ)i + A + M , (2) 



where 1 is the unit matrix, A is antisymmetric and M is symmetric and traceless. With (|^) equation (|^) reads: 
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where the first term is the isotropic Heisenberg interaction, the second one - the Dzyaloshinskii-Moriya interaction and 
the last one - the anisotropic pseudo-dipole interaction. The vector Di^ contains the three independent components 
of Aij. In general the occurrence of Dij requires low crystal symmetries. The origin of the anisotropic interaction 
Di j and Mij is the spin-orbit coupling. In the last few years there is some effort concerning the question of including 
the term with A/^.j in the calculations |^^, but in what follows we neglect Mij, because it is only of second order in 
the spin-orbit coupling constant A whereas Dij is of first order in A. 

The weak ferromagnetic moment in the predominantly antiferromagnetic ordered CuO planes in La2Cu04 can be 
interpreted with a small spin canting. This spin canting could be described via a DM interaction term. In a 
general form of the DM vector is introduced and it is shown, that only a vector Di j which varies from bond to bond 
corresponds to the crystal structure and is able to describe the observed weak ferromagnetism. In (2jj2^ Coffey and 
coworkers consider different crystal symmetries and get Di j for the whole lattice by requiring for the vector Di j that 
the energy of any configuration of spins is invariant under the symmetry transformations of the crystal structure. 
Once Di j is fixed on a single bond, the symmetries determine Di j on the entire lattice. 

Following 1^,^ we consider in this paper two different symmetries for Dij which correspond to the orthorhombic 
and tetragonal phases of La2Cu04 ||30|-p^. The arrangements of the atoms in a CuO-plane are presented in figure 1. 
For the orthorhombic arrangement (figure la) one finds the following Di^: 

• Da,b = (di,d2,0) 

• Da,c = {-(h-, -di,0) 

• Db.D = ~Da,C = (^2,^1,0) 

• Dc,D = -Da,b = (-di,-d2,0) 

• Dc,E = -Da,c = (rf2,rfl,0) 



For the tetragonal arrangement (figure lb) the Dij is as follows: 

. Da^b = (0,di,0) 
. Da^c = (0,d2,0) 

• Db,d = ~Da,c = (0,-^2,0) 

• Dc,D = ~Da,b = (0,-di,0) 
. Dc,E = -Da,c = (0,-^2,0) 

Furthermore we will distinguish between equal and different signs of the parameters di and ^2 . Following the arguments 
of we will restrict our considerations to di and d2 of equal strength. Let us define four different cases for the 
parameters: orthorhombic symmetry with either di ~ +d2 or di = —d2 and tetragonal symmetry with either di = +d2 
or di ~ —d2. 

First let us consider the classical Hamiltonian 

N 

-ffm/ = ^ A.j(5,®5j) (4) 

with the configurations for Dij as given above. Si and Sj are classical 3D vectors and the summation runs over 
nearest neighbour bonds. The classical groundstate of (^) was discussed in detail in ||2^. Here we briefly summarize 
some ground state features which are relevant for the further discussion. The energy of any nearest neighbour bond 
(i-ij) on the lattice can be minimized by spin vectors Si and Sj perpendicular to each other and perpendicular to the 
DM vector Dij; the corresponding optimum bond energy is Eij — — \ Dij \ S*^. Because of the special symmetry 
I di 1 = 1 d2 I there is no frustration in the Hamiltonian (|^) and the total groundstate configuration can be build by 
the suitable arrangement of optimized bonds. In table 1 we present some important groundstate features for the 
considered four different cases of Dij. 
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case 


crystal DM 
structure vector 






net spin 
FM configuration 


la 


orthorhombic di — 


-d2 


= d 


no 


4 SL 


lb 


orthorhombic di — - 


-d2 


^ d 


yes 


2 SL 


Ila 


tetragonal di = 4 


-d2 


= d 


yes 


2 SL 


lib 


tetragonal di = - 


d2 


= d 


no 


4 SL 



Table 1 : The four considered cases with some of its ground- 
state features of the classical DM Hamiltonian Hdm 

Here FM means ferromagnetic moment and 2 SL or 4 SL means 2 or 4 sublattices, respectively. As an example the 
groundstates for the cases la and lb are illustrated in more detail in figure 2. 
In the next section we turn to the quantum system. 

III. THE MODEL 

We start with the so called Ji — J2 model on the square lattice: 

N 

1=1 

Si denotes the spin-i operator on site i and x and y are the unit lattice vectors in x- and y-direction. 

An intensive investigation over the last years (see e.g. 0- suggests the following phase diagram: For small 
values of J2/J1 the system shows AFM-LRO. At some critical value J2"* / Ji of about ~ 0.4 the conventional coUinear 
AFM-LRO breaks down. Then a region of a spin-liquid state could be realized for 0.4 < J2/J1 < 0.65. For values of 
J2/ Ji > 0.65 an AFM-LRO arises within the initial sublattices (4-sublattice AFM). 

Now we add the Hamiltonian of the DM-interaction 

H ~ Hj^j^ + Hdm (6) 



TV 



(7) 



(8) 



Using spin-flip operators we can rewrite the DM term 

Dij{Si X Sj) — —\^iD^j [S'j Sj — S^Sj — S^Sj + S^S^] 

+ ^Dt^[StS- ^SIS+]] 

The Hamiltonian Hdm contains terms of the form S^S^ which do not commute with the z-component of the total 



spin S. Hence we have to take all 2^-Ising states for the construction of the wave function. We use a modified 
Lanczos procedure to calculate the groundstate of the system. 



IV. NUMERICAL RESULTS 



Based on the properties of the pure Ji — J2 model we discuss several order parameters, cf. which describe 

the relevant magnetic properties of the model. 

For dominating Ji the corresponding antiferromagnetic LRO parameters are the components of the square of the 
staggered magnetization: 



1 ^ 



7 = 



1 for ieA 
— 1 for ieB 



(9) 
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where is the 7-component of the spin on site i, Ti is the corresponding staggered factor and A and B are the two 
sublattices. Obvious these parameters describe the ordinary 2-sublattice antiferromagnetic ordering. 
For dominating J2 the relevant order parameters are 



^ N/2 



^^x,y,z ; Ti^Q = ±1 ; a = A, B . (10) 



The sum runs over the sites of the sublattice a. These order parameters describe antiferromagnetic ordering within 
the two initial sublattices A and B, when the whole system forms a 4-sublattice AFM. For the rotational invariant 
model without DM interaction the x,y,z components of the respective parameters are identical, but the DM interaction 
breaks this rotational symmetry. 

For D ^ a weak ferromagnetism may occur. The relevant parameters for describing this weak ferromagnetism 
are the components of the square of the magnetization, 



1 ^ V' 

i=l 



7 = . (11) 



Obviously (Af)^ is the expectation value of the square of the total spin (^S^ 

In the pure model {Dij — 0) the Hamiltonian (^) commutes with and and the groundstate is isotropic with 
= (and consequently (A/^)^ — 0). If we now put on the anisotropic DM interaction {Dij ^ 0) the Hamiltonian 
does not commute with 5^ and we may expect (Af')^ 7^ 0. Koshibae et al pointed out that anisotropic spin 
correlations are a result of the interplay between quantum fluctuations and the anisotropic DM interaction. This can 
be illustrated as follows: Consider a classical spin system without DM interaction and a coUinear spin structure. The 
DM term contains the cross product of the spins so it vanishes. If now the quantum fluctuations are put on, the DM 
interaction favours the direction of spin fluctuations with minimal energy. Of course this direction depends on the 
DM vector but additionally on the classical coUinear structure. 

We have calculated the groundstate for a square lattice with = 16 and A^ = 20 sites for the four different 
configurations of the DM vector given in table 1. In the following we focus our interest on the orthorhombic symmetry. 



since this symmetry is valid for low doping in La2-2:Ba2,Cu04 32 1. Furthermore we present preferably the data for 
A^ ~ 20. The corresponding data for A^ = 16 are qualitatively the same. For simplicity we put Ji = 1 for the rest of 
the paper. 

First we consider the situations of J2 = and of J2 = 1. Then the isotropic model {Dij = 0) is coUinear long-range 
ordered with either a 2-sublattice structure (J2 = 0) or a 4-sublattice structure (J2 — 1). The influence of the DM 
interaction on the magnetic order parameters is shown in figures 3 and 4. 

Figure 3 shows the different components of the ferromagnetic moment (Af^)^ as a function of the DM interaction 
for J2 = 0. The results are as follows: The DM interaction causes a small ferromagnetic moment {{M'^)'^ ^ 0) which 
increases with d. This moment is anisotropic: ((Af^)^ > (A^^'^)^. The ferromagnetic moment is maximal for case 
lb which is in accordance with the classical situation (see table 1). We have calculated the same parameters for 
J2 = 1. In that case the classical structure is a 4-sublattice AFM, which does not have a net ferromagnetic moment. 
Consequently we find in the quantum case that the ferromagnetic moment is about 10 times smaller as for J2 = 0. 

In the next two figures 4a and 4b we present the different components of the AFM-LRO parameters. For J2 = the 2- 
sublattice AFM-LRO parameters {M])^ have to be considered and for J2 = 1 the 4-sublattice AFM-LRO parameters 
(-^^a) ^'"^ relevant. As discussed above the anisotropy depends on the DM vector Dij and on the underlying 
magnetic structure. Correspondingly the anisotropy is changed going from J2 = to J2 = 1 or, alternatively, going 
from case la to case lb. Furthermore, in dependence on the symmetry of the DM vector, we can observe either an 
enhancement or a suppression of a certain component of the antiferromagnetic order parameter. For instance the 
z-component of the 4-sublattice AFM-LRO parameter {M^ ^ is enhanced by the DM interaction if Dij favours 

a classical 4-sublattice structure (see figure 4a) whereas the same order parameter is suppressed if Dij favours a 
classical 2-sublattice structure (see figure 4b). 

The comparison of the magnitude of ferromagnetic and antiferromagnetic order parameters (figures 3 and 4) shows 
that the ferromagnetic moment remains small even if the DM interactions reaches the same strength as the exchange 
coupling. 

Let us now discuss the antiferromagnetic order parameters in dependence on the frustration parameter J2/J1. For 
this we select from figure 4 in each case the dominating components (Af^)^ and (Af^^)^. The magnetic order for 
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Di j = is determined by the competition between Ji and J2 . For dominating Ji and for dominating J2 the magnetic 
structure is cohinear antiferromagnetic with two or four sublattices. The transition between both phases is in the 
classical case at J2 — ^ Ji and may be connected with a spin liquid phase in the quantum case. Since the DM term 
also favours a two (case lb) or a four sublattice structure (case la) we expect that this competition is influenced if we 
put on the DM term. The result should be a shift of the transition either to lower or higher values of J2. 

As shown in figure 5 we find indeed the expected shift, it is in particular strong for case la. As discussed above, 
the critical J2 for the pure model {Dij=0) is about J2 ~ 0.4, which is too large to be realistic for the cuprate 
superconductors. We conclude that a DM interaction of this symmetry is able to support the breakdown of the 
antiferromagnetic long-range order due to frustration and to shift the critical J2 to more realistic values. (Notice, that 
the abruptness of the change in the order parameters seen in figure 5 as well as in figure 6 is due to the well-known 
level crossing for the iV = 20 lattice (see e.g. g, [0)). 

Finally we will discuss the region of a possible spin liquid in more detail. As pointed out in there are two 

interesting candidates for unconventional noncoUincar ordering in this region. One is a vector chiral order parameter 
introduced in 0,0, 



/1 r 1 ^ 

\ L i=l 

Cli = 8e^ ((5, X S,) + {S, X Si) + {Si x 5,)) . 



(12) 
(13) 



This order parameter measures the handiness of a plaquette of three spins. For the isotropic Heisenberg system this 
parameter is the same for all three components x,y and z. But in our case of anisotropic DM interaction the different 
components may have different values. 

Another candidate for exotic ordering is the dimer (or spin-Peierls) order parameter. 



D 



1 ^ 

-y(-i) 

1=1 



' SiSi-]-g 



(14) 



which measures the long-range phase coherence of spin dimers (singlets of two spins) along the x (or y) direction. 

In figure 6 we present these parameters C and D, where we have selected the largest component of C. Both 
order parameters show in the pure model a characteristic maximum at about J2 = 0.55. The influence of the DM 
interaction is twofold. First the maxima are shifted in the same way as discussed for figure 5. Second the maxima are 
suppressed particularly for D. The suppression of the maximum in D can be simply understood by having in mind 
that the dimer ordering is connected with a singlet groundstate which is obviously not realized for Di j ^ 0. On the 
other hand, the effect of the DM interaction on the chiral order parameter seems to be more complicated. The 
DM term creates a spin canting which might under certain circumstances support the realization of a vector chiral 
ordering. Indeed, we found one example {N = 16, di = +d2 — 0.3), where the chiral order parameter is enhanced 
in a small region around J2 ~ 0.4 (figure 6c). Here the DM interaction term together with the frustrating J2 causes 
some canted groundstate with enhanced vector chirality. However, it remains unclear whether this increase of is 
an artifact of the particular cluster symmetry of N = 16. 



V. SUMMARY 



The magnetic properties of the Ji — J2 model with anisotropic Dzyaloshinskii-Moriya (DM) interaction are de- 
termined by the competition between terms favouring a coUinear 2-sublattice AFM and terms favouring a coUinear 
4-sublattice AFM. The former one is stabilized by the Ji term and the latter by the J2 term. The DM term may favour 
the 2-sublattice structure as well as the 4-sublattice structure in dependence on the symmetry of the DM vector Dij 
(see table 1 in section 2). As a result the DM interaction shifts the spin liquid region (which separates the two coUinear 
antiferromagnetic phases for small and large J2) to larger or smaller values of J2. Because the critical J|''** w 0.4 for 
the breakdown of the 2-sublattice LRO in the pure model {Di,j=0) is much larger than realistic J2 values for cuprate 
superconductors fl-^ additional mechanisms are needed to shift the transition to smaller values of J2 . Recently it 
has been shown that static holes simulating doping can shift J|''** in the desired direction. In this paper we find 
that a DM term with suitable symmetry acts in the same direction: In case of orthorhombic symmetry a DM vector 
with identical x- and y-components (case la in table 1) is appropriate to weaken the 2-sublattice antiferromagnet and 
to realize lower J|''**. 
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Besides of the shift of the transition region the DM term creates an anisotropy in the spin correlations. Which 
components are enhanced and which are suppressed is a resuh of the interplay between quantum fluctuations, the 
particular symmetry of the DM vector as well as the underlying magnetic structure of the pure model. Hence it occurs 
that for the 2-sublattice AFM another component is favoured than for the 4-sublattice AFM. 

The third remarkable effect of the DM interaction is the appearance of a weak ferromagnetic moment. For the 
considered finite lattices this moment occurs for all symmetries of Dij. In accordance with the classical picture the 
orthorhombic structure with different x and y components of Di j (case lb in table 1) yields the strongest ferromagnetic 
moment and might indicate the existence of weak ferromagnetism in the thermodynamic limit just for this symmetry 
of Dij. 

Finally we analyzed exotic dimer and vector chiral order parameters showing a characteristic maximum in the spin 
liquid region. This maximum is in general suppressed by the DM term. 

Besides of the most important orthorhombic symmetry considered preferably in the paper we have also calculated 
the magnetic quantities for the tetragonal symmetry (cf. figure lb and table 1). In principle the same scenario as for 
the orthorhombic symmetry is valid, i.e. in dependence on the strength and the symmetry of the DM vector Dij we 
may have anisotropy in the spin correlations, a net ferromagnetic moment and a shift of the spin liquid phase. 
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FIG. 1. The orthorhombic (a) and tetragonal (b) arrangements of the CuO-octahedra taken from [24]. The filled circles are 
the copper sites which carry the spins, the open circles are oxygen sites, which are tilted up out of the plane and the crossed 
circles are oxygen sites, which are tilted down out of the plane. The square represents one tilted CuO-octahedra. The points 
A-E denote certain copper sites used in the text. The arrows show one particular arrangement of the DM vector Di,j for this 
lattice (cases la and Ila, see table 1. 

FIG. 2. The classical spin configurations for the Hamiltonian (|^; (a): case la with a 4-sublattice configuration and without 
a net ferromagnetic moment, (b): case lb with a 2-sublattice configuration and with a net ferromagnetic moment. (All spins 
are aligned either along the z-axis or in the xy-plane.) 

FIG. 3. The x,y and z components of square of the magnetization (M^)^ versus d/Ji for case la (full line) and case lb 
(dashed hue) for iV = 20 and J2 = 0. 

FIG. 4. The 2-sublattice AFM-LRO parameters {M^f for J2 = (fuU hues) and the 4-sublattice AFM-LRO parameters 
{M2,a) for J2 = 1 (dashed lines) versus d/Ji for A*' = 20; (a): case la, (b): case lb . 

FIG. 5. The dominant 2-sublattice AFM-LRO parameters {M^f (fuU lines) and the dominant 4-sublattice AFM-LRO 
parameters (MJ^q) (dashed lines) versus J2/J1 for N = 20 and different values of d; (a): case la with (M^'^f and (M,%) , 
(b): case lb with {M^f and (A^f;i')^ 

FIG. 6. Exotic order parameters versus J2/J1 for different d: (a) - The dominant chiral order parameter C^'^ for case la 
(solid line) and for case lb (dashed line) for A*' — 20; (b) - The dimer order parameter D for case la (solid line) and case lb 
(dashed line) for A = 20 ; (c) - The dominant chiral order parameter C^'^ for case la (solid line) and C for case lb (dashed 
line) for N = 16. 
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